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CONSIDERATION OF SIMPLE MODEL ATOMS FOR 


ARGON DISCHARGE CALCULATIONS* 
by James R. Rose and John V. Dugan, Jr. 

Lewis Research Center 

SUMMARY 

Simple model atoms are useful for satisfactory calculation of plasma properties in 
gas discharges with low electron temperature T g . Model atoms of 2-, 3-, and 5-levels 
give ionization fractions f in satisfactory agreement with 26-level results for cesium 
plasmas that are optically thick to resonance radiation. The calculation of nonequili- 
brium f values is extended to argon discharges of T g values of 8000 and 11 605 K. 

The calculated values of f closely approximate the Saha predictions in the completely 
thick radiation limit. Level inversions are calculated at both T g values for argon 
plasmas that are optically thick to resonance radiation. 


INTRODUCTION 

Numerical values for free-electron distribution functions f e (u) and electronic state 
populations have been obtained for cesium-seeded argon discharges with low elec- 
tron temperature T g from simultaneous solution of the Boltzmann and rate equations 
(ref. 1). (Symbols are defined in the appendix. ) The rate equations were solved for a 
5-level cesium (Cs) model atom in reference 1, the ground electronic level plus four 
excited states. The limitations of 2-, 3-, and 5-level model atoms in predicting de- 
partures of electron number densities N g from equilibrium (Saha) values were reported 
in reference 2. These results for an assumed Maxwellian f e (u) were compared with 
results for a 26-level model (ref. 3). Where comparison is possible, the agreement is 
within a factor of 2 at low and high initial cesium number densities Nq s but is only 

* Presented at the Second Arc Symposium sponsored by the American Physical 
Society, Hartford, Conn. , October 20, 1970. 


22 ? 

qualitative at intermediate and low cesium nuclei densities N^ g (<10 m ) for 
plasmas optically thick to resonance radiation. 

The electron distribution functions calculated for cesium in reference 1 are near 
Maxwellian down to T g = 2100 K for plasmas optically thick to resonance radiation. 
Numerical difficulties were encountered for T g < 2100 K, where the iterative scheme 
for f e (u) does not converge for the high-energy tail. 

Simple model atoms serve two useful purposes: they allow for (1) convenient cal- 
culation of nonequilibrium N g values in low T g plasmas and (2) straightforward isola- 
tion of numerical difficulties that may occur in the solution for the troublesome high- 
energy tail of f e (u) (refs. 1 and 2). 

The purpose of this work is to investigate the usefulness of relatively simple (2-, 

3-, and 5 -level) argon model atoms in calculating N g values for argon discharges. 

The electron distribution function is assumed to be Maxwellian. The ratio of the first 
excited state to ground state binding energy of argon is much different from that of 
cesium (0.26 and 0.64, respectively). This essential difference in the atomic structure 
of argon and the current interest in argon discharges warrants this separate investiga- 
tion. The investigation entails calculating N g values in argon discharges operating 
with electron temperatures of 8000 and 11 605 K (1 eV) for the 2-, 3-, and 5-level 
model atoms. The calculations were done for plasmas that are (1) thin to all radiation, 
(2) thick to resonance radiation, and (3) thick to all radiation. The thick to all radiation 
limit should approximate the Saha values of N g . 

ATOMIC MODELS 

The 2-, 3-, and 5-level model atoms are shown schematically in figure 1. The 5- 
level model consists of a ground state and four excited states. The lowest lying excited 

5 

state is obtained by grouping all of the energy levels arising from the 3p 4s electronic 
configuration (i.e., the 4s p^ and 4s pg ^ energy levels). Similarly, the second 
state is obtained by grouping the levels arising from the 3p^4p electronic configuration 
and the third by grouping those states arising from the 3p^5s and 3p®3d configurations. 
Finally, all the remaining states are approximated by considering a continuum of states 
from the highest lying level of the 3p 5 3d state to the ionization limit. The energy of 
the fifth level is taken to be the mean energy of this upper continuum of states. The 
degeneracy of the fifth level was chosen in the same manner as in the cesium studies 
(ref. 2). This level must approximate the presence of several thousand spectroscopic 
states in the absence of comparison with results for more detailed models. The Saha 
check calculations (completely thick to all radiation) and the nonequilibrium results 
of this report are quite insensitive to the degeneracy of the fifth level gg. Values of 
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50, 200, and 500 were tried for gg. The equation N g = 0 was best satisfied for 

gg = 200, so the results are reported for this value. This is consistent with the manner 

in which the cesium results of reference 2 varied with the upper-level degeneracy. 

The energy of the upper level of the 3-level atom is taken to be the mean energy of 

5 

the highest lying level of the 3p 4p electronic configuration and the ionization limit of 
argon. Its degeneracy is taken to be that of the 3p^4p electronic configuration. The 
2-level atom consists simply of the ground state and the 3p^4s states of argon. 

The binding energy of a given level is assigned as the mean energy of the multiplet 
it represents, that is, 


E 


L 




(la) 


and the corresponding degeneracy is simply 


g L 2 g Lj 
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Values of the E^j and g-^j 


were obtained for argon from reference 4. 
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ASSUMPTIONS AND ANALYSIS 

The equations solved in the present work are the steady- state atomic rate equations 
and the mass conservation equation. The rate equation for the Lth state, that is, the 
rate of change of the number of bound electrons in level L, N^, is (ref. 1) 

N l = - N L (N e< r L ♦ + N e N k K k _ l 

K*L 

+ y. »A.r. * fe aP T + h 1 (m‘ 3 ) (sec- 1 ) (2) 

K>L 
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In the steady state N L = 0. 

In equation (2) K K _ L ((m 3 ) (sec”'*')) is the rate coefficient for populating the Lth 
state by collisional excitation or deexcitation of the Kth state; 


M T 


.(■!> 

\ K*L 


v- . Tfionl 

jL— K + K L 
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((m ) (sec )) is the rate coefficient for depopulating the Lth 


-L 


state by collisional excitation, deexcitation, and ionization; Aj^jJsec” ) is the proba- 
bility of a spontaneous radiative transition from the Kth state to the Lth state (K>L); 

■E Al— K/ * s *** e probability of spontaneous radiative transitions from the Lth 
K<L / * 

state to lower states K; (K^ ap ) ((m 3 ) (sec”*)) is an effective two-body rate coefficient 
for capture into level L; and /3^ ((m 3 ) (sec”*)) is the rate coefficient for radiative cap- 
ture into level L. Photoionization is ignored because of its unlikelihood. The mass 
conservation equation can be written approximately as 



N o ■= N i + Z N L <”" 3 > < 3 > 

L 

where N Q is the initial argon number density. Equations (2) and (3) are solved for the 
electronic state populations and the electron number density N g (= N^) for assigned 
values of both the argon initial number density N Q and the electron temperature T g . 

Boundary and spatial effects are neglected in the calculations and it is assumed that 
the electron distribution function is Maxwellian at the assigned temperature T . The 
later versions of the Gryzinski inelastic cross sections (ref. 5) are used to calculate the 
excitation and ionization rate coefficients. The corresponding inverse processes are 
then calculated by the principle of detailed balancing. 


Radiative Transition Probabilities 

The resonance transition probability A 2 ~l was obtained from reference 6. 
Approximate, scaled values of were used for the other resonance transitions, 

that is, 


a k-l = 


g 2 <E 1 - E K> . 

; 2 2-i 

g K (Ej - E 2 ) 2 


(sec *) 


(4) 
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for K > 2 . To obtain this expression it was assumed that the oscillator strengths for 
all resonance transitions are nearly the same and equal to the oscillator strength for 
the lowest resonance transition. 

The value for A^g is the mean value of the thirty allowed transitions between 
the 3p^4p and 3p^4s multiplets (ref. 7). Other transition probabilities were either 
assigned approximately from comparison with similar transitions or calculated using 
the Coulomb approximation (ref. 8). A summary of the numerical values used for the 
A k _^ l is given in table I. 


Radiative Capture Coefficients 

The radiative capture coefficients for the excited states (L > 1) were calcu- 
lated by detailed balancing using hydrogenic photoionization cross sections (ref. 9). 

The ground state value of the radiative capture rate coefficient was obtained 
from an experimental ground state photoionization cross section (ref. 10). This photo- 
ionization cross section was curve fitted against incident photon energy, and /3^ was 
obtained from it by using the principle of detailed balancing. 

The variation of the ionization fraction f in cesium plasmas has been discussed in 
terms of the source and sink terms in the electron conservation equation (ref. 3). The 
normalized expression for N 0 = 0 is 

M .MM 

f2 X! ( K L P ) = Z *£% (5) 

L=1 L=1 L=1 

where f = N e /N 0 is the fraction ionized and n^ = N l /N q is the normalized population 

of the Lth electronic state. The quantity K^ n ((m 3 ) (sec - *)) is the rate coefficient 
for ionization from the Lth state. Although the equation for the ionization fraction, 
equation (5), appears to be simply quadratic in f, the normalized bound state populations 
n-^ are strong functions of N Q and optical thickness for a given electron temperature. 

The levels which are sources of ionization for the various model atoms can be de- 
termined from the net ionization rates. The net ionization rate from the Lth electronic 
state can be written from equation (5) as 

R net = K^ n n L - |(k° ap )* f 2 + /3 L fj (m 3 ) (sec -1 ) (6) 

n pt 

The positive Rl values define states which are free electron sources and the negative 
net •*-* 

values define states which are free electron sinks. The behavior of the fraction 
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ionized as 


in 


T>net 

*1, * 


T and N are changed can generally be explained by observing the changes 

w U 


RESULTS AND DISCUSSION 
Saha Check in Limit of Complete Optical Thickness 

It is possible to perform a numerical check of the accuracy of a model atom in the 
limit of complete optical thickness. In this limit the true fraction ionized is that pre- 
dicted by the Saha equation for argon, including all possible occupied electronic states. 
The results of this check showed that the 2-, 3-, and 5-level model atoms predict the 
Saha result to within 5 percent for all cases tried (T = 8000 K, 11 605 K, and 
N q = 10 18 to 10 26 ). 


Optically Thin 

The calculated fraction ionized for an optically thin plasma and a 1-eV electron 

temperature is plotted against initial argon number density in figure 2. The figure 

shows that the fraction ionized is generally much less than the Saha solution, although 

the 5-level model atom results approach the Saha results at the highest N values. The 

18 22 8 ® 

fraction ionized is constant for N Q values of 10 to 10 m for the 5-, 3-, and 2- 

level model atoms. The 5-level model atom predicts a sharp increase in the fraction 

22 -3 

ionized as N Q increases above 10 m . The 2- and 3 -level model atom predictions 
of the fraction ionized remain relatively constant over the entire range of N Q values 
considered. 

A better understanding of these fraction- ionized variations is obtained from exam- 

«pt 

ining figure 3. The figure shows the positive R T values given by equation (6) scaled 

•*-* Mpf 

by a quantity R milv , which is the maximum sum obtained by adding the positive Rv 

1 1 iii ax 

values for each N value for the 5-level model atom calculations for given plasma 

o n i. 

radiation conditions. The quantity R^ /R max * s denoted by Rjy This scaling by 

R „„ correlates the behavior of the 5-level model atom results to the 2- and 3-level 
max 

results. States that are free electron sinks are not shown in the figure. The apparent 
discontinuity in some of the curves is due to the lack of data between integral powers of 
N Q ; that is, a given curve starts at the lowest integral power of N Q for which it has a 
positive value. 

As can be seen in figure 3, the value of RI is roughly constant and equal for all 

18 A 22 

three model atoms in the N Q range of 10 to 10 , where the fraction ionized behaves 
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similarly (fig. 2). A look at the rates involved shows that in this range of N Q the frac- 
tion ionized predicted by the three model atoms considered corresponds to the coronal 
limit (ref. 11), that is, 


f = 


M 

E 

L=1 

M 

Ea 

L=1 


(7) 


It can be seen from equation (7) that generally in the coronal limit radiative capture 
balances collisional ionization. In the case of argon at these conditions only the ground 
state participates, so that 


f 


K 10n n 
K 1 n l 


h 


( 8 ) 


The reason for the sharp increase in f for the 5-level atom at higher N Q values is 

apparent from figure 3. All states except the ground state contribute increasingly to the 

free-electron population as N Q increases until maxima in the values (2 < L < 5, 

M = 5) are reached at the same value of N that f has a maximum. Also note the 

o 04 -3 

upper level is the main ionization source at the higher N Q values. For N q > 10 m 

three-body capture becomes a more dominant process which causes all levels except 

the L = 5 level to become electron sinks. It can be seen in figure 2 that there is a 

corresponding decrease in the fraction ionized. 

It is clear that the upper level is populated by a ladder mechanism (ref. 2) in the 
5-level case. Also, this fifth state has a relatively low binding energy (0. 78 eV), so 
that once it is significantly populated it is readily ionized since kT g > Eg for a 1-eV 
electron temperature. For 2- and 3-level model atoms there is no sharp increase in 
the fraction ionized as N Q increases, and the upper level never becomes a source of 
free electrons. The energy gap between the second and third levels of the 3-level atoms 
is so large that the collisional excitation rate of the third state from the second state is 
much less than the depopulation rate of the second state by spontaneous transitions to 
the ground state A 2 _^. Thus, the ladder mechanism is inhibited for the simpler model 
atoms. 

The results for T 0 = 8000 K, optically thin to all radiation (not shown) are roughly 
the same for the 5-, 3-, and 2-level atoms. The fraction ionized is roughly constant 
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over the entire range of N Q values considered (f = (5 to 8)xl0 ). The f curves re- 

main flat because the excitation rate coefficients for populating even the lowest lying 
excited states are small compared to the corresponding values at kT g = 1 eV. Thus, 
the states which are relatively easily ionized once populated never become significantly 
populated. 


Optically Thick to Resonance Radiation 

Population inversion . - For plasmas optically thick to resonance radiation the phe- 
nomenon of population inversion (n^ > n^) occurs for both electron temperatures con- 
sidered (8000 and 11 605 K). The normalized excited state number densities n^ for 
the 5-level atom are displayed against N in figure 4 for kT = 1 eV (11 605 K) and in 
figure 5 for T g = 8000 K. The figures also show the curve of ^ against N Q calcu- 
lated for thick to all radiation conditions ng Saha' can seen i n figures 4 and 5, 
the level inversion (n^ > n^) occurs in the range where ng is greater than or not much 
less than ng Saha’ ^ hen the population of the L = 2 level is this large, it can readily 
populate states above it by excitation. Excitation from the L = 2 state is the largest 
contributor to the populations of both the L = 3 and L = 4 states. The third level 
(L = 3), however, is more readily depopulated than the L = 4 level by spontaneous 
radiative transitions. The other depopulating mechanisms for these states are relatively 
small at the lower N Q values. Thus, as the L = 4 and L = 3 states become popu- 
lated, the L = 3 state is preferentially depopulated by spontaneous transitions, and 
this depopulation in turn leads to n^ > n^ over a range of N Q values. 

As N q increases, however, the electron density N g (= fN Q ) increases, depopula- 
tion by spontaneous transitions is overshadowed by collisional depopulating mechanisms, 
and the population inversion no longer occurs. Since the collisional depopulating mech- 
anisms depend on N g , one might expect the population inversion to occur up through 
higher N Q values for T g = 8000 K (where the fraction ionized is relatively small) than 

for kT^ = 1 eV (11 605 K). Figures 4 and 5 show this to be the case. 

e 23 -3 

Figure 4 shows n^ > n^ for N Q values greater than 10 m . This evident level 

inversion is probably due to the artificial manner in which the L = 5 level is con- 
structed, and there is no physical interpretation given for this. 

Fraction-ionized variation . - Energy equivalent of electron temperature, 1 eV 
(11 605 K): The fraction ionized is plotted against initial argon number density for 
plasmas thick to resonance radiation in figure 6 for kT g = 1 eV (11 605 K). The f 
values predicted by the 2- and 5-level atoms are roughly the same at low and high N Q 

values. At intermediate N values the 5-level atom predicts f values up to 50 percent 

° 24 -3 

higher than either the 2- or 3-level atom. At the higher N values (>10 m ) three- 
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body capture becomes the dominant capture process which causes the fraction ionized to 

decrease. This decrease in f is much more pronounced for the 3-level atom than for 

either the 2- or 5-level atom. All model atoms predict fractions ionized much below 

the Saha values except at the highest N Q values where the dominant atomic processes 

are collisional and Saha equilibrium is approached. 

The f curves for the 2- and 3-level atoms can be predicted by equation (7) (thick 

18 22 -3 

coronal limit) in the range N Q = 10 - 10 m . In particular, though, 


f 


_ «-ion 
n 2 K 2 

h 


in this range of N Q . 

Figure 7 shows the normalized ionization rates against the initial argon number 
density for the various model atoms for a 1-eV electron temperature and plasmas thick 
to resonance radiation. These correspond to the fraction-ionized curves shown in fig- 
ure 6. In the case of the 5-level atom all states except the ground state are contributors 
to the free-electron population over most of the range of N Q values, and the normalized 
ionization rates roughly follow the trends of the fraction-ionized curve (fig. 6). Similar 
remarks apply to the normalized ionization rate curve for the 2-level model atoms. 

It is seen from figure 7 that the ionization rate curve for the 3 -level atom increases 
as N q increases for N Q > 10 m . This obviously does not explain the rapid de- 
crease in the fraction ionized for the 3-level atom in the same region (fig. 6). The 
reason for this rapid decrease in f is that in this region the dominant atomic rate is 
the three-body capture rate into the upper level of the atom. 

Electron temperature, 8000 K: The fraction ionized is plotted against initial argon 
number density in figure 8 for T g = 8000 K and the plasmas optically thick to resonance 
radiation. Qualitatively the results are similar to the corresponding results for kT g = 

1 eV (11 605 K) except that the curves are shifted to higher N Q values. 

The normalized ionization rates plotted against initial argon number density for the 
8000 K electron temperature (fig. 9) show the same shift relative to the corresponding 
kT = 1 eV (11 605 K) results (see fig. 7). 


Effect of Optical Transparency 

Figures 10 and 11 show the effect of the optical transparency of the plasma on the 
fraction ionized. The curves are the ratios of the calculated values of the fraction 
ionized for the 5-level atom to the fraction ionized predicted by the Saha equation for 
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electron temperatures of 8000 and 11 605 K, respectively. The curves for plasmas 
thick to all radiation are very nearly horizontal and show the accuracy of the 5 -level 
model atom in this completely thick limit. Comparison of figures 10 and 11 shows that 
for a given plasma optical thickness the plasma is farther from equilibrium at the lower 
temperature. The figures also indicate a trend toward equilibrium as N q increases, 
as should be the case. 


CONCLUDING REMARKS 

Calculations made using 2-, 3-, and 5- level argon model atoms indicate that these 
simple model atoms are useful in predicting the fraction ionized in argon discharges. 
All three model atoms predict the expected coronal limit at low initial argon number 
densities N Q . At high N Q the results for the 5-level atom approach the Saha limit, 
which is the limiting solution for this condition. The accuracy of these simple model 
atoms at intermediate N Q values is difficult to assess because of the lack of other data 
for comparison. 

The 5-level argon model atom proved useful in interpreting level inversions that 
occurred at both temperatures used in the calculations (8000 and 11 605 K) for plasmas 
thick to resonance radiation. 

Simple model atoms are also useful for isolating numerical difficulties that may 
occur in the numerical solution of the Boltzmann equation in regions where inelastic 
collision effects are not negligible. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 29, 1971, 

129-02. 
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APPENDIX - SYMBOLS 


probability of spontaneous radiative transition L— K, sec 

•^L 2 A L-K 

K<L 



binding energy of level L, eV 
fraction ionized 

electron energy distribution function 
degeneracy of level L 

3 -1 

rate coefficient for collision-induced transition L— K, (m ) (sec ) 

3 - 1 

effective two-body rate coefficient for capture into level L, (m ) (sec ) 

3 - 1 

rate coefficient for ionization from level L, (m ) (sec ) 

I K L-K+4° n 


K*L 

kT g energy equivalent of electron temperature, eV 
N number density 

n^ dimensionless number density of Lth level 

R™ ax maximum ionization rate from level L for 5-level model atom, (m^) (sec 

net ionization (capture) rate of level L, (m^) (sec - *') 

Rj^ scaled net ionization rate, R^ e ^/R™ ax 

T g electron temperature, K 

/3-j^ rate coefficient for radiative capture into level L, (m ) (sec ) 

Subscripts: 

Cs cesium 

e electron 

i ion 

K, L atomic level indices 

M upper level 

max maximum 


mu 


i 


o initial number density 
Superscripts: 
cap capture 
ion ionization 

time derivative 
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Figure 7. - Normalized net ionization rate from Lth 
state as function of initial argon number density. 
Energy equivalent of electron temperature, 1 eV 
(11 605 K); plasma thick to resonance radiation. 
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Initial argon number density, Nq, m" 3 

Figure 8. - Fraction ionized as function of initial argon 
number density. Electron temperature, 8000 K; plasma 
thick to resonance radiation. 



Figure 9. - Normalized net ionization rate from Lth state 
as function of initial argon number density. Electron 
temperature, 8000 K; plasma thick to resonance radia- 
tion. 



Ratio of fraction ionized to Saha value, f/^saha 



Figure 10. - Ratio of calculated fraction ionized 
to Saha fraction ionized for various plasma 
optical conditions. Electron temperature, 
8000 K. 



Figure 11. - Ratio of calculated fraction ionized 
to Saha fraction ionized for various plasma 
optical conditions. Energy equivalent of elec- 
tron temperature, 1 eV (11 605 K). 
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